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Numerical  study  of  a  loss-cone  negative  mass  instability  to  amplify  whistler  waves  by  energetic 
electrons  in  the  radiation  belts  is  presented.  The  results  show  that  a  very  low  intensity  whistler  wave 
can  be  amplified  by  50  keV  electrons  more  than  25  dB,  consistent  with  the  Siple  experimental  result 
[Helliwell  et  al. ,  J.  Geophys.  Res.  85,  3360  (1980)].  The  dependencies  of  the  amplification  factor 
on  the  energetic  electron  density  and  on  the  initial  wave  intensity  are  evaluated.  It  is  shown  that  the 
amplification  factor  decreases  as  the  initial  wave  intensity  increases.  However,  this  gain  can  still 
exceed  15  dB  for  a  30  dB  increase  of  the  initial  wave  intensity,  which  is  needed  for  the  purpose  of 
precipitating  MeV  electrons  in  the  radiation  belts.  We  then  show  that  there  exists  a  double  resonance 
situation,  by  which,  as  an  example,  a  wave  is  simultaneously  in  cyclotron  resonance  with  50  keV 
electrons  as  well  as  with  1.5  MeV  electrons;  the  wave  is  first  amplified  by  50  keV  electrons  and  then 
precipitates  1 .5  MeV  electrons.  With  the  aid  of  the  cyclotron  resonance,  the  threshold  field  for  the 
commencement  of  chaos  in  the  electron  trajectories  is  reduced  considerably  from  that  for  a  general 
case.  Pitch  angle  scattering  of  1.5  MeV  electrons  is  demonstrated.  The  results  show  that  a  whistler 
wave  with  magnetic  field  amplitude  of  0.08%  of  the  background  magnetic  field  can  scatter  electrons 
from  an  initial  pitch  angle  of  86.5°  to  a  pitch  angle  <50°.  ©  2007  American  Institute  of  Physics. 

[DOI:  10.1063/1.2743618] 


I.  INTRODUCTION 

In  the  magnetosphere,  energetic  electrons  in  the  radia¬ 
tion  belts  are  trapped  by  the  Earth’s  dipole  magnetic  field 
and  undergo  a  bouncing  motion  about  the  geomagnetic  equa¬ 
tor.  Very  energetic  electrons  (those  in  the  MeV  levels)  have  a 
strong  impact  on  passing  satellite  systems.  Satellites  are  de¬ 
signed  to  survive  a  certain  amount  of  radiation  (ionizing) 
dose  accumulated  during  their  lifetimes.  Unexpected  en¬ 
hancement  of  the  radiation  fluxes  caused  by,  for  example, 
very  strong  solar  storms,  will  significantly  increase  the  total 
radiation  dose  to  the  satellites.  Consequently,  the  radiation 
damage  on  active  electronics  and  detectors  of  satellite  sys¬ 
tems  will  accumulate  much  faster  than  designed  for.  As  the 
damage  exceeds  a  threshold  level,  satellite  systems  become 
incapable  of  performing  their  mission.  Therefore,  mitigation 
of  unexpected  radiation  enhancement  is  an  issue  of  great 
concern. 

Propagation  of  whistler  waves  in  the  magnetosphere  has 
been  studied  extensively.1  Such  waves  can  be  excited  during 
lightning  events2,3  or  launched  into  the  magnetosphere  by  the 
ground4  or  space-based  VLF  transmitters.  Whistler  waves 
can  be  ducted  in  an  L-shell  of  the  magnetosphere  to  continu¬ 
ously  interact  with  the  energetic  electrons  trapped  in  the 
same  L-shell.  The  motions  of  very  energetic  electrons  are 
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adversely  affected  by  the  wave  fields  which,  in  essence,  scat¬ 
ter  some  that  are  trapped  in  the  radiation  belts  into  a  loss 
cone.1  Induced  electron  precipitation3  5-8  by  whistler  waves 
has  been  observed. 

The  Doppler  shifted  electron  cyclotron  resonance 
interaction9-14  has  been  suggested  to  be  responsible  for  the 
correlation  of  the  presence  of  magnetospheric  whistler  waves 
and  the  electron  precipitation  events.  However,  the  simulta¬ 
neous  observation  of  precipitation  events  at  geomagnetic 
conjugate  regions  in  both  hemispheres  due  to  a  single  light¬ 
ning  flash3  requires  a  different  theoretical  interpretation. 
While  the  numerical  results14  show  that  the  electron  cyclo¬ 
tron  resonance  interaction  is  indeed  effective  in  scattering 
very  energetic  electrons  into  the  loss  cone,  in  particular, 
those  electrons  with  their  initial  pitch  angles  close  to  the  loss 
cone  angle,  the  percentage  of  the  total  number  of  very  ener¬ 
getic  electrons  which  can  be  resonant  with  the  wave  at  a 
given  frequency  is  small.  It  suggests  that  this  process  (rely¬ 
ing  on  cyclotron  resonance  interaction)  would  require  that 
the  wave  have  a  broad  frequency  spectrum,  so  that  a  substan¬ 
tial  fraction  of  the  very  energetic  electrons  can  be 
precipitated. 

In  our  early  work,15-17  we  have  shown  that  the  trajecto¬ 
ries  of  trapped  very  energetic  electrons  in  the  presence  of  a 
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whistler  wave  can  become  chaotic.  Once  it  occurs,  many  of 
them  will  wander  into  the  loss  cone.  This  chaotic  scattering 
process  can  be  an  effective  approach  for  controlling  the 
population  of  very  energetic  electrons  in  the  radiation  belts. 
There  is  a  threshold  field  requirement  for  the  commencement 
of  chaos  in  the  electron  trajectories.  Favorably,  it  has  been 
observed18  20  that  trapped  energetic  (keV)  electrons  in  the 
magnetosphere  can  significantly  amplify  whistler  waves.  An¬ 
other  supporting  evidence  of  whistler  wave  amplification  by 
energetic  electrons  is  the  natural  event  of  chorus21,22  occur¬ 
ring  in  the  inner  magnetosphere.  The  energy  transfer  is  also 
through  the  electron  cyclotron  resonance  interaction.23-26 
This  suggests  an  optimal  approach,  which  applies  the  chaotic 
scattering  process  under  a  double  resonance  situation,  for  the 
control  of  the  population  of  very  energetic  (MeV)  electrons 
trapped  in  the  magnetosphere.  This  approach  is  first  to  am¬ 
plify  the  incident  whistler  waves  by  the  trapped  keV  elec¬ 
trons  (having  a  loss-cone  velocity  distribution)  in  the  radia¬ 
tion  belts;  the  amplified  whistler  waves  then  scatter  the 
trapped  MeV  electrons,  which  are  also  in  cyclotron  reso¬ 
nance  with  the  wave,  into  the  loss  cone.  In  the  previous 
work,27  5  keV  electrons  were  considered  for  wave  amplifi¬ 
cation  and  double  resonance  conditions  could  not  be  satisfied 
simultaneously.  Thus,  the  chaotic  scattering  process  did  not 
invoke  cyclotron  resonance;  consequently,  its  threshold  wave 
field  was  shown  to  be  quite  high. 

In  the  present  work,  we  show  that  when  50  keV  elec¬ 
trons  are  considered  for  the  wave  amplification,  wave  can 
also  be  in  cyclotron  resonance  with  1.5  MeV  electrons  dur¬ 
ing  the  pitch  angle  scattering.  Thus  the  threshold  for  the 
commencement  of  chaos  is  reduced  considerably.  In  Sec.  II, 
this  amplification  process  is  analyzed  numerically.  Amplifi¬ 
cation  gain  vs  the  energetic  electron  density  is  evaluated.  The 
peak  intensity,  the  modulation  period,  and  the  gain  of  the 
amplified  wave  vs  the  incident  wave  intensity  are  also  deter¬ 
mined.  In  Sec.  Ill,  pitch  angle  scattering  of  very  energetic 
electrons  by  a  whistler  wave  under  double  resonances  is 
studied.  We  first  use  the  double  resonance  conditions  to  de¬ 
termine  the  initial  parameters  of  those  very  energetic 
(1.5  MeV)  electrons  which  are  resonantly  interacting  with 
the  wave.  These  parameters  are  then  used  in  the  evaluation 
of  pitch  angle  scattering.  Summary  and  conclusions  are  pre¬ 
sented  in  Sec.  IV. 

II.  LOSS-CONE  NEGATIVE  MASS  INSTABILITY: 
AMPLIFICATION  OF  WHISTLER  WAVES 

The  magnetospheric  electron  plasma  consists  of  three 
components:  (1)  cold  background,  (2)  energetic  (in  the  keV 
range)  electrons,  and  (3)  very  energetic  (in  the  MeV  range) 
electrons.  The  background  electron  plasma  is  nonrelativistic 
(y~  1)  and  the  momentum  P  of  the  electron  has  a  loss  cone 
distribution,  which  can  be  modeled  as 

Fb(P{  ,PZ)  =  nb(7T-mlj\)(lP]b)-{1^]pV 

Xexpf-  (I*  +  P$l\P)bl  (1) 

where  j  is  the  loss  cone  index,  and  the  loss  cone  angle  ipL 
=  tan “ 1 V7 ;  &Pjb=[mTeb/{\/2+j/3)]m\  subscript  b  denotes 
“background.” 


The  background  electron  plasma  determines  the  propa¬ 
gation  characteristics  of  the  whistler  wave,  which  has  the 
dispersion  relation 

(o=n0c2k2/w2pb,  (2) 

where  copb=(4Trnbe2/m)m\  nb  is  the  background  cold  elec¬ 
tron  density;  il0=eB0m/c  is  the  nonrelativistic  electron  cy¬ 
clotron  frequency. 

Energetic  electrons  carry  weak  relativistic  effects  (y0 
— 1.1,  i.e.,  50  keV  electrons)  and  have  significant  density  nE 
to  amplify  whistler  waves.  The  electron’s  momentum  also 
has  a  loss  cone  distribution  given  by 

fE(PlyPz)  =  27rP_FE(P^Pz) 

=  ne(2iTin/jl)(\Pjt)-aj*3)P2{*' 

Xexp[-  (P\  +  P2Z)I  \P2r] 

=  nJel(PJfcz(Pz),  (3) 

where  \Pje=[mTet/(\/2+j/3)]in;  Ttr>Teb,  subscript  e. 
stands  for  “energetic.” 

In  the  following,  an  instability  process  to  amplify  whis¬ 
tler  waves  by  energetic  electrons  having  a  loss  cone  velocity 
distribution  in  the  magnetosphere  is  studied.  The  loss  cone 
distribution  makes  this  component  potentially  unstable  to  mi¬ 
croinstability.  The  relativistic  effect,  through  the  Doppler 
shifted  cyclotron  resonance  interaction  provides  an  effective 
channel  for  the  energy  exchange  between  the  electrons  and 
the  waves.  The  resonance  condition  is  given  by 

(o=  ildy+  k  •  v,  (4) 

where  a><  n0  for  the  whistler  waves;  o>  and  il0  are  the  wave 
frequency  and  the  nonrelativistic  electron  cyclotron  fre¬ 
quency,  and  y  is  the  relativistic  factor  of  the  electron.  For 
small  y,  i.e.,  o><fl 0/y,  the  resonant  electrons  are  moving 
oppositely  to  the  wave  propagation  direction. 

The  relativistic  effect  provides  nonlinearity  in  wave- 
electron  interaction,  which  results  in  bunching  those  elec¬ 
trons  in  near  cyclotron  resonant  interaction  with  the  wave. 
The  electron  bunching  provides  a  positive  feedback  to  excite 
the  loss-cone  negative  mass  instability  for  the  wave 
amplification.28  In  the  present  work  we  consider  only  wave 
amplification  in  time,  a  local  uniform  magnetic  field  will  be 
assumed.  The  resonant  trajectory  equations  of  a  single  elec¬ 
tron  in  the  whistler  wave  fields  and  in  a  uniform  background 
magnetic  field  z  B0  are  given  by 

da/dr  =  -  (k/il0)(eE0/m)(  1  +  kvjio) cos  </>,  (5) 

d</>/d/  =  -  Aw+  (k/il0)(eEo/m)(\  -k-kvJ(o)a~l  sin  </>,  (6) 

dvJdt  =  ( (Iq/ k) {eEJmc2) ( <*/  y2) (kc2/ a>  +  iOcos  </>,  (7) 

dy/df  =  -  (flo/^)(^£,0/mc2)(tt/y)cos  </>,  (8) 

where  a=kv j/Hq  and  <£=&Zo+[0o+<p(ri]--/oA<i/dt' ;  A (o 
=  \a)0+{l0{y-y0)/yy0+k(vz-vz0 ),  A  ojq  =<*>-<«*),  and  <**> 
=  Il0/  yo-kv^,  00=tan~l(v>^/vxo)i  ip  accounts  for  the  phase 
shift  in  the  electron  gyration  due  to  interaction  with  the  wave 
fields,  y0  is  the  initial  relativistic  factor  of  the  resonant  elec- 
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FIG.  1 .  Momentum  distribution  of  energetic  electrons  and  the  regions  close 
to  cyclotron  resonant  interaction  with  the  wave. 


trons.  The  ratio  of  (7)  and  (8)  leads  to  an  invariant  relation 
y(kc2/ u)+vz)= const  =  y0(kc2/ (o+v#),  which  is  used  to  ob¬ 
tain  &to=&<oQ-(k2c2/ (o-coQ)(y-y0)/ y.  It  is  noted  that 
k2c2 / uj- d)0> 0  for  whistler  waves  considered  in  the  present 
work,  i.e.,  Aa>  increases  as  electrons  lose  energy  to  the  wave 
and  vice  versa. 

Velocity  of  the  electron  plasma  in  the  magnetosphere  has 
a  loss-cone  distribution.  Only  a  small  portion  of  the  total 
number  of  electrons,  e.g.,  in  the  two  shade  regions  in  the 
velocity  distribution  shown  in  Fig.  1,  is  near  the  resonant 
interaction  with  the  wave.  The  wave  is  experiencing  cyclo¬ 
tron  damping  to  the  electrons  which  are  initially  at  exact 
cyclotron  resonance  with  the  wave,  i.e.,  Aa>o=0.  However, 
the  wave  can  exchange  energy  with  other  electrons  having  a 
mismatch  frequency  Aojq  slightly  different  from  zero.  De¬ 
pending  on  the  initial  phases  of  those  electrons,  the  interac¬ 
tion  can  cause  them  either  to  gain  energy  from,  or  lose  en¬ 
ergy  to,  the  wave,  initially. 

In  region  I,  Aoj01<0;  electrons  reduce  their  mismatch 
frequencies  in  the  interaction  as  losing  energy  to  the  wave 
(y-  y0<0);  it  results  in  the  increase  of  the  interaction  period 
of  losing  energy  to  the  wave.  In  contrast,  as  electrons  gain 
energy  from  the  wave  (y-yo>0),  their  mismatch  frequen¬ 
cies  increase  and  interaction  periods  of  gaining  energy  from 
the  wave  decrease.  On  average,  those  electrons  with  Aco^ 
<0  lose  energy  to  the  wave.  For  those  electrons  in  region  II 
with  A<o0=A<oo2>0,  the  above  dynamic  interaction  process 
is  reversed,  and  on  average,  those  electrons  gain  energy  from 
the  wave.  Using  the  definition  Ao)0=(o-aj0=(o-(}0/ y0 
+  kvz o,  Aa>oi<A(o02  leads  to  uzoi<uzo2’  which  indicates  that 
there  are  more  electrons  in  the  Ao>01<0  region  than  in  the 
Ao)02>0  region  (i.e.,  N{>N2  as  shown  in  Fig.  1).  Overall, 
the  wave  gains  energy  from  electrons  and  thus  it  is  amplified. 

The  amplification  is  the  result  of  the  collective  effect  in 
the  electron-wave  resonance  interaction.  Thus  the  single  par¬ 
ticle  resonant  trajectory  equations  (5)-(8j  are  averaged  over 
the  electron’s  random  phase  angle  (with  respect  to  the  wave 
field)  to  obtain  the  governing  equations  for  the  induced  cur¬ 
rent  densities  in  the  energetic  electron  plasma.  The  phase 
average  is  based  on  a  similar  procedure  as  that  presented  in 
Ref.  27.  After  phase  average,  (8)  is  converted  to  an  energy 
conservation  equation,  which  can  be  integrated  to  obtain  the 
relation  ( A  o>)K ,  =  A  <doi  +  0)[eo(  1  +  er)  /2  y0AwKmc2] 

x[£q(/)-£q(0)],  where  ( )e|  represents  an  average  over  the 
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initial  random  phases  of  those  energetic  electrons  in  region  I 
of  Fig.  1;  -N2  (**-  !)/«(***.) 

-/£l(P  Li)/gZ(Pz,)](2yimA^//:)3  is  the  net  density  of  ener¬ 
getic  electrons  transferring  energy  to  the  wave  through  the 
resonant  interaction  with  the  wave;  er=\+co2pbl  id{\q 
=  (kc/o))2  is  the  dielectric  function  of  background  plasma 
responding  to  the  whistler  wave. 

The  phase  averaged  equations  derived  from  the  three 
first  order  differential  equations  (5)-(7)  are  then  combined 
with  the  wave  equation,  which  is  a  second  order  differential 
equation,  to  yield  a  single  differential-integral  equation,27  in 
essence  a  fifth  order  ordinary  differential  equation 

[d]  +  2Yd]  +  (Aw2  +  C)d,]E0 

=  a0{l+(Nc/\nE)(io2pl/a0^l) 

x[e0(l  +  Er)IAyaLnFmc2\El(t)  -  E2(0)]} 

x  f  E0(/')cos<A<#/-/'))e,d/\  (9) 

Jo 

where  r  =  -[3fir/8(l  +Er)](bxaJ  (x>)  and 

Aa)2  =  Acoq,  +  «q{52[(1  +  Vfld >/kc2)  +  3k(co/c)2(\  +  kvz q/cj) 

-  (m)a2k(u>/c)4(\  +  kc^co)2]E2 

-  a>Aa>o( l  +  fcu^o/w)} 

-  (3fc,  «o/4^  -  2^sE2)AoIsE0- 

s  =  (k/il0)(e/m)(  1  +  tojo/wJrtQ1 , 

/c=  f  E0(t' -  f'))Eidt' 

Jo 

and 

/s=  f  E0(Osin(A<Mf-f’)>,.|dt'; 

Jo 

C=  [Aa)2f/(1  +er)y0]((l  +  kv^oj)  -  \  12(^1^  y^cc)2], 

a0  =  -  [Ao)p(  n2A<o0la^/(  1  +  e.r)  yow], 
and 

A  (o2e  =  A  nt:e2/me0. 

Equation  (9)  governs  the  temporal  evolution  of  the  field 
amplitude  of  the  whistler  wave  interacting  resonantly  with 
this  particular  group  of  trapped  electrons  (in  the  shade  re¬ 
gions  of  Fig.  1)  in  the  radiation  belts.  Equation  (9)  is  linear¬ 
ized  to  obtain  the  relation  |Aojo,|3  =  /^(o2  filial/ 16(  1 
+  er)  yo<*>,  which  is  used  to  determine  |Ao>qi|.  We  next  analyze 
this  equation  numerically  to  determine  the  conditions  for  the 
amplification  of  whistler  waves. 

For  numerical  analysis,  (9)  is  normalized  to  a  dimen¬ 
sionless  form  by  introducing 

X=[e0(fir-  l)|Atuoi|/4a)y0Awe;nc2]l/2Eo(rt. 
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FIG.  2.  Temporal  evolution  of  the  amplitude  of  a  whistler  wave  traversing 
the  magnetosphere. 


f  =  l^oik* 
and 

g  =  (NJ\nt)(er+  \)kc(ow2pbHl0(er-  l)A<Uo, 

=  (Nt/bne){er+l)(a>/Lu>0l)2. 

The  normalized  equation  has  the  form 

[d\  -  48  Xfj]  +  ( 1  +  bX2  +  48 fJOdfVC  =  1 6[  1  -  g(X2  -  X2)]fc 

(10) 

where  X0=X(0);  /c=/£X(f  )cos  )d£'  and  fs 

=-J(oX((')  sin  rff.f  W ;  <p(t,£')=ff,{l-2g[X2(£") 

-xftW; 

b  =  \4\o)2prk2c2/y0U:r-  l)|AoJoi|3<»ir(l  +  v^to/kc2) 

+  3(a>IVro*2c2)  -  ( 9/4 )  «g(  ro*2c2  )2] 

=  192[(er+  1  )/(er  —  l)](wro/«o“o); 

b  and  g  are  constant  coefficients;  (10)  is  subjected  to  the 
initial  conditions:  X(0)=Xo,  dfX( 0)=v3X0,  d2eX( 0)  =  3X0. 

In  the  present  work,  we  consider  the  resonant  energetic 
electrons  of  50  keV,  which  is  different  from  the  previous 
case27  for  5  keV.  This  significant  increase  of  the  energy  of 
the  resonant  energetic  electrons  is  necessary  in  order  to 
achieve  the  double  resonances  as  being  shown  in  the  next 
section.  Having  the  background  parameters,  b- 1440  and  g 
=  2X  107  (i.e.,  er~4  and  IAcoqJ/w— 5  X  10"4),  and  setting 
X0=3.58  X  10~4,  (10)  is  solved  by  an  ordinary  differential 
equation  (ODE)  solver.  The  result  is  presented  in  Fig.  2, 
showing  the  temporal  evolution  of  the  field  amplitude  E0(t). 
The  dB  scale  plot  in  Fig.  2  is  for  a  direct  comparison  with 
the  early  Siple  experimental  result.18-20  In  Siple 
experiments,18  injected  Siple  signals  of  3  kHz,  propagating 
along  the  7,  =  4  shell,  were  often  amplified  by  10-30  dB  and 
oscillated  in  time.  The  numerical  result  presented  in  Fig.  2 
also  shows  that  whistler  wave  amplitude  is  amplified  more 
than  25  dB  by  trapped  relativistic  electrons  through  the  loss- 
cone  negative  mass  instability  and  also  oscillates  in  time  in  a 
similar  fashion,  a  good  agreement  with  the  experimental 
observation. 

The  dependence  of  the  wave  amplitude  gain  G  on  the 
density  of  energetic  electrons  is  through  the  parameter  g  in 


-  Log  (g) 


FIG.  3.  Amplitude  gain  G  of  the  amplified  whistler  wave  as  a  function  of 
the  density  of  energetic  electrons  (*-log(g)). 


(10).  Since  g  is  inversely  proportional  to  (A ne/Ne) 
X  (Aci>oi/a>)2  *  An^/3/A^e,  the  gain  G  in  dB  is  plotted  against 
-log(g)  in  Fig.  3,  where  G=20  log(£0m/^in)  and  E0m  and  Ein 
are  the  maximum  amplitude  of  the  amplified  wave  and  the 
amplitude  of  the  incident  wave,  respectively.  As  shown,  the 
gain  G  increases  (nearly  linearly)  with  the  density  of  ener¬ 
getic  electrons,  as  expected. 

It  is  noted  that  the  field  amplitude  Ein  of  the  incident 
whistler  wave  in  Fig.  2  is  rather  low  (to  be  consistent  with 
that  of  previous  experiments);  however,  in  the  practical  ap¬ 
plication  for  achieving  significant  electron  precipitation,  the 
incident  wave  field  has  to  increase  considerably.  Therefore,  it 
is  important  to  realize  how  the  gain  G  varies  with  the  inci¬ 
dent  wave  intensity  for  a  fixed  background  condition.  This  is 
exemplified  by  considering  a  case  with  the  following  back¬ 
ground  conditions:  O0/a>=2,  y0=l.l,  and  er=ll;  lAa^J/io 
~5X  10”4  and  u/c=0.417;  setting  v±/c= 0.32,  leads  to  a0 
=0.358,  6=456,  and  g  =  3.2xl06.  The  peak  intensity 
7m(*£o/J  °f  the  amplified  wave  versus  the  intensity  /in 
(*£2n)  of  the  incident  wave  is  presented  in  Fig.  4.  As  shown, 
the  rate  of  increase  of  Im  decreases  as  /in  increases.  In  other 
words,  G  (7in)  is  a  decreasing  function  of  7in,  which  is  pre¬ 
sented  in  Fig.  5.  It  shows  that  with  a  30  dB  increase  of  7in,  G 
is  reduced  by  about  10  dB.  The  process  of  wave  amplifica¬ 
tion  is  via  phase  bunching  of  gyrating  electrons  as  the  posi¬ 
tive  feedback  of  wave-electron  interaction.  However,  when 
the  phase  bunching  exceeds  180°,  the  interaction  turns  to 
debunch  electrons.  Thus  the  amplitude  of  the  wave  oscillates 
during  the  continuous  interaction.  The  oscillation  period  is 
the  sum  of  the  bunching  and  debunching  times,  which  de¬ 
crease  as  the  wave  intensity  increases.  Consequently,  the 
modulation  period  T  (in  seconds)  of  the  wave  intensity  is 
expected  to  decrease  with  7in  because  7m  increases  with  7in. 
This  dependence  T  (7in)  is  also  presented  in  Fig.  4. 
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FIG.  4.  The  dependency  of  the  peak  intensity  Im  and  the  modulation  period 
T  of  an  amplified  wave  on  the  intensity  7jn  of  the  incident  wave. 


III.  DOUBLE  CYCLOTRON  RESONANCES 
FOR  EFFECTIVE  PRECIPITATION  OF  VERY 
ENERGETIC  ELECTRONS  BY  WHISTLER  WAVES 

Whistler  waves  can  cause  the  trajectories  of  trapped  en¬ 
ergetic  electrons  in  the  radiation  belts  to  be  chaotic,  leading 
to  a  scattering  of  very  energetic  (MeV)  electrons  into  the  loss 
cone.  However,  this  chaotic  scattering  process  requires  a 
threshold  wave  field  for  the  commencement  of  chaos  in  the 
electron  trajectories.  The  amplification  process  considered  in 
the  prior  section  is  aimed  at  reducing  the  power  requirement 
of  the  incident  whistler  waves.  In  addition,  cyclotron  reso¬ 
nance  is  found  to  be  the  optimal  condition  to  achieve  the 
minimum  threshold  power  requirement.  To  take  the  advan¬ 
tage  of  both  effects,  a  double  resonance  condition  has  to  be 
satisfied,  namely,  the  wave  can  be  simultaneously  in  cyclo¬ 
tron  resonance  with  energetic  (a  few  tens  of  keV  level)  elec¬ 
trons  as  well  as  with  very  energetic  electrons. 

A  wave  resonant  with  MeV  electrons  satisfies  a  reso¬ 
nance  condition,  similar  to  that  given  by  (4), 


(o  =ilo/y2  +  kvZ2-  (11) 

To  show  that  a  double  resonance  situation  is  possible,  we 
first  express  the  resonance  condition  in  a  general  form  w 
=  n0ly+kPz/ym,  where  y=(\+P\lm2c2+P2zlm2c2)m ,  P 
=  ymv L,  and  Pz=ymvz.  Because  of  the  y  dependence,  this 
condition  leads  to  a  quadratic  equation  for  Pz  as  AP2 
+  2BPZ+C=0 ,  where  A  =  (l -o>2/*2c2),  £=miy*,  and  C 
=  {m/k)2[£ll-(i)2(\  +P2  /m2c2)].  This  quadratic  equation  has 
two  real  solutions  Pz=[-B±{B2-AC)in]/Ay  subject  to  the 
condition  B2^AC.  The  double  solutions  suggest  that  the 
wave  can  be  simultaneously  resonant  with  two  different 
groups  of  electrons.  The  coefficients  A  and  C  of  the  quadratic 
equation  are  positive  because  co/kc<  1  and  i\!  yx%2>  co  are 
considered;  thus  both  Pz  are  negative,  i.e.,  the  two  groups  of 
electrons,  which  can  resonantly  interact  with  the  wave,  move 
opposite  to  the  propagation  direction  of  the  wave. 

We  now  find  the  initial  conditions  of  electrons  such  that 
(4)  and  (11)  can  be  satisfied  simultaneously.  The  relation¬ 
ships  and  notations  to  be  applied  are  first  introduced  as  fol¬ 
lows:  n0*£"3;  poj6*L”3/2;  mb*  (1  -  y^Ym£lQlL  is  the 
bounce  frequency  of  trapped  very  energetic  electrons;  n 
=  (pu)2bl  o>fl0)1/2  is  the  index  of  refraction  of  the  background 
cold  plasma  on  the  whistler  wave;  L  is  the  number  of  earth 
radii,  i.e.,  L  value  of  a  magnetic  flux  tube;  0,=tan_I(/),/2  is 
the  initial  pitch  angle  of  the  energetic  electrons  which  con¬ 
tribute  to  wave  amplification,  thus  cos  0,=  1  /(I  +y)1/2;  yx 
~  1.1  is  the  relativistic  factor  of  the  energetic  electrons;  02  is 
the  initial  pitch  angle  of  trapped  very  energetic  electrons 
which  are  being  precipitated;  the  normalized  phase  velocity 
of  the  wave  V-vpl c=(fl0/ (i>p)£~xl2,  where  (=il0/(o. 

With  the  aid  of  some  known  background  parameters: 
nb~280  cm"3  at  L= 4.9,  0.25  Gauss  at  L=l,  and 

ayn0~(l/3)X  10"2  for  y2=3  and  L= 2,  we  can  obtain  the 
functional  dependence  of  the  background  parameters  on  L  as 
P<o*=2ttX  1.63  X  106/L3/2,  n0=2TrX7X  ltf/L3,  and 
<oBnx=(u)B/u))i(ii0/a))=i.oix  io-3(i-y22)“1/2/£. 

The  resonance  conditions  (4)  and  (11)  and  the  dispersion 
relation  lead  to 


no/<0=  72  +  ,*(y2-  1 ) 1/2  COS 

=  yx  +  /i(y[-  1) 1/2 cos  0j,  (12) 

n  =  2.33L3/2£l/2  =  2.33L3/2[y2  +  /i(y^  -  1)1/2  cos  02}m . 

(13) 

Equations  (12)  and  (13)  are  solved  to  yield 

cos  02=  l/2[(yj-  l)/(?2-  1)]1/2  X  {{y2lyx  +  l)cos  0, 

“  (yjyx  -  l)[cos2  01  +4y,/5.43L3(yf-  1)]1/2}. 


(14) 


FIG.  5.  Amplitude  gain  G  of  the  amplified  whistler  wave  as  a  function  of 

the  intensity  of  the  incident  whistler  wave.  Thus,  (12)  and  (13)  become 
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Hq/cj  ->  n0  =  [  y2( yf  -  1  )m  cos  0\ 

-  7i(?2-  1) 1/2  cos  l)l/2  COS  0, 

-(>2-  1) 1/2  cos  02}n 
=  (?2-  7iVt(7|-  l)l/2cos  0, 

-(■)!-  1) 1/2  cos  02 ] 

and  the  normalized  (to  cu)  bounce  frequency  of  the  trapped 
very  energetic  electrons  and  the  normalized  (to  c)  phase  ve¬ 
locity  of  the  wave  are  given  by  a>fl=7.07  X  10"3(1 
-y2 2r1/2n0/L  and  V=0.432L"3/2ao1/2. 

The  electron  trajectory  equations  are  similar  to  the  set  of 
equations  used  to  derive  the  resonant  trajectory  equations 
(5)— (8 ),  except  an  equivalent  (restoring)  force  to  simulate  the 
magnetic  mirror  effect  is  also  included  in  the  equations.  After 
a  canonical  transformation  to  a  new  set  of  canonical  coordi¬ 
nates,  the  normalized  electron  trajectory  equations  for  the 
normalized  canonical  coordinates  (P,Q,z,pz)  in  the  whistler 
wave  fields  are  given  by17 

dP/dt  =  -  y2  -  l)G  +  (nofli/y2)sinz,  (15) 

dQ/dt=(il0/y2 -  l)(P/a0)  +  ayy2)cosz,  (16) 

dz/dt  =  p:/y2,  (17) 

dpjdt  =  -  o)2bz  +  (fly/y2)(P  sin  z  +  il0Q  cos  z),  (18) 

where  y2 = k(  a >/c)[k(c/ cd)2 + P2 + p2z  +  Ogg2  +  ft  2 

+  2ft,(P cos  z-ft0G  sin z)]1/2;  £lx=eB{lm  and  B ,  is  the  am¬ 
plitude  of  the  wave  magnetic  field.  In  terms  of  the  canonical 
coordinates,  the  new  Hamiltonian  is  found  to  be  time  inde¬ 
pendent,  i.e.,  a  constant  of  motion.  In  essence,  (15)-(18)  de¬ 
scribe  trajectories  in  a  three-dimensional  phase  space.  It  is 
also  noted  that  if  coB  is  set  to  zero,  (15),  (16),  and  (18)  can  be 
combined  to  obtain  an  additional  constant  of  motion,  \{P2 
+  ft202)-ftopz=const.  This  constant  reduces  the  system  to 
two  autonomous  equations  and  the  trajectory  will  have  regu¬ 
lar  behavior  in  a  two-dimensional  phase  space.  In  other 
words,  the  bouncing  motion  (i.e.,  the  magnetic  mirror  effect) 
is  essential  in  inducing  chaos  in  the  electron  trajectory.  This 
is  because  bounce  causes  the  phase  relationship  between  the 
wave  and  the  electron  to  change  abruptly,  so  the  electron 
trajectory  can  evolve  from  a  regular  behavior  to  a  chaotic 
one.  However,  the  chaoticity  of  the  system  is  not  very  sen¬ 
sitive  to  the  value  of  the  bounce  frequency,  as  long  as  it  is 
not  too  small.17 

We  now  consider  the  y2=4  case  that  precipitates 
1.5  MeV  electrons.  From  (14),  02=86.5°.  We  then  have  ft0 
=  8.6038,  V=0.0517,  and  a^=2.945  X  10"2.  Equations 

( 1 5)— ( 1 8)  are  now  integrated  numerically  to  evaluate  the 
pitch  angle  scattering,  resulting  from  the  wave-electron  reso¬ 
nance  interaction.  Presented  in  Fig.  6  is  a  result  that  double 
resonance  condition  is  satisfied.  Thus  the  result  presented  in 
Fig.  6  is  different  from  those  presented  in  Fig.  5  of  Ref.  16 
and  in  Fig.  4  of  Ref.  17.  As  shown,  when  the  wave  magnetic 
field  B i  increases  to  7X(1.16X  10“4  B0 ),  large  pitch  angle 
scattering  occurs.  In  other  words,  with  wave  magnetic  field 
amplitude  of  0.007/8.6038=0.0008=0.08%  of  the  back- 
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FIG.  6.  (Color  online)  Temporal  evolution  of  the  pitch  angle  of  a  1.5  MeV 
electron,  which  is  interacting  with  a  whistler  wave  at  Doppler  shifted  cyclo¬ 
tron  resonance.  The  wave  magnetic  field  increases  consecutively  from  1 
to  7  X  (1.16  X  10"4  B0). 


ground  magnetic  field,  electron  is  scattered  to  a  pitch  angle 
<50°,  less  than  the  loss  cone  angle.  Further  increase  of  the 
wave  magnetic  field  B{  to  exceed  7X(1.16X  10~4  B0 ),  i.e., 
H,  >0.007,  the  numerical  results  confirm  that  large  pitch- 
angle  scattering  persists.  This  is  an  example  that  a  whistler 
wave,  with  proper  parameters,  can  be  amplified  by  the  elec¬ 
trons  having  energies  in  tens  of  keV  ranges  and  simulta¬ 
neously  scatters  MeV  electrons,  both  processes  via  the  cy¬ 
clotron  resonance  interaction.  Comparing  the  result 
presented  in  Fig.  6  with  those  of  the  nonresonant  cases  pre¬ 
sented  in  Refs.  16,  17,  and  27,  the  cyclotron  resonance  inter¬ 
action  reduces  the  threshold  wave  magnetic  field  amplitude 
for  achieving  precipitation  by  more  than  an  order  of  magni¬ 
tude.  However,  the  required  interaction  time  is  increased  by 
more  than  two  orders  of  magnitude. 

IV.  SUMMARY  AND  CONCLUSION 

It  is  shown  that  energetic  electrons  in  the  radiation  belts 
can  amplify  a  whistler  wave  through  a  loss-cone  negative 
mass  instability.  The  gain  factor  increases  nearly  linearly 
with  the  density  of  energetic  electrons.  The  results  show  that 
this  amplification  process  can  enhance  the  field  intensity  of  a 
very  low  amplitude  incident  whistler  wave  by  more  than 
25  dB.  The  theory  is  formulated  and  the  numerical  result  is 
shown  to  agree  well  with  the  experimental  result,  both  quali¬ 
tatively  and  quantitatively.  However,  the  peak  amplitude  of 
the  amplified  wave  still  depends  on  the  amplitude  of  the 
incident  wave.  In  order  to  achieve  significant  electron  pre¬ 
cipitation,  the  field  amplitude  has  to  exceed  a  threshold  level. 
Hence,  the  amplitude  of  the  incident  wave  has  to  increase 
considerably  from  the  level  of  the  previous  Siple  experi¬ 
ments,  as  well  as  from  those  of  currently  available  VLF 
transmitters.  It  is  found  that  this  gain  decreases  as  the  inci¬ 
dent  wave  intensity  increases;  however,  it  reduces  to  less 
than  10  dB  for  a  30  dB  increase  of  the  incident  wave  inten¬ 
sity.  In  other  words,  this  amplification  process  will  still  be 
effective  and  reduce  considerably  the  required  field  intensity 
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of  the  incident  whistler  wave  for  the  purpose  of  precipitating 
those  very  energetic  electrons  in  the  MeV  range. 

An  optimal  approach  to  reduce  the  population  of  MeV 
electrons  trapped  in  the  magnetosphere  is  suggested  as  fol¬ 
lows.  To  precipitate  those  electrons  with  pitch  angles  close  to 
the  loss  cone  angle,  the  required  wave  field  is  below  a  thresh¬ 
old  for  the  onset  of  chaos.14  Small  pitch  angle  scatterings  via 
electron  cyclotron  resonance  interaction  will  be  adequate  to 
diffuse  those  electrons  into  the  loss  cone.13  However,  to  pre¬ 
cipitate  those  electrons  trapped  far  away  from  the  loss  cone, 
it  will  need  a  chaotic  scattering  process.  In  this  case  the  wave 
field  has  to  exceed  a  threshold  for  being  able  to  cause  the 
electron  trajectories  chaotic.  Although  the  threshold  field  is 
reduced  considerably  at  cyclotron  resonance,  the  numerical 
result  shows  that  it  is  still  quite  large.  The  numerical  results 
demonstrate  that  a  1 .5  MeV  electron  can  be  scattered  from 
an  initial  pitch  angle  of  86.5°  to  a  pitch  angle  <50°  by  a 
whistler  wave  with  the  magnetic  field  amplitude  of  0.08%  of 
the  background  magnetic  field.  It  converts  to  about  3100  pT 
at  L= 2,  and  200  pT  at  L=5.  Thus  wave  amplification  is  in¬ 
deed  needed.  We  have  shown  that  the  background  energetic 
electrons  (e.g.,  <100keV  electrons)  can  amplify  injected 
whistler  waves  through  a  loss-cone  negative  mass  instability 
by  more  than  15  dB.  We  have  further  shown  that  a  double 
resonance  condition  can  be  satisfied,  namely,  that  the  wave 
can  be  cyclotron  resonant  with  energetic  (a  few  tens  of  keV 
level)  electrons  for  amplification,  and  yet  still  be  able  to 
precipitate  MeV  electrons  via  effective  cyclotron  resonance 
interaction.  Finally,  it  should  be  pointed  out  that  this  optimal 
approach,  relying  on  electron  cyclotron  resonance  interac¬ 
tion,  requires  that  the  wave  have  a  broad  frequency  spec¬ 
trum,  so  that  a  considerable  fraction  of  the  very  energetic 
electrons  can  be  precipitated  simultaneously.  A  broadband 
wave  can  be  considered  as  a  sum  of  many  narrow  band  com¬ 
ponents  and  different  components  have  different  central  fre¬ 
quencies.  Each  narrow  band  wave  component  will  interact 
with  a  distinctive  group  of  electrons  determined  by  the  cy¬ 
clotron  resonance  conditions.  Thus  the  amplification  and  pre¬ 
cipitation  processes  for  each  narrow  band  wave  component 
will  be  similar  to  those  for  a  monochromatic  wave  consid¬ 
ered  in  the  present  work.  However,  when  the  wave  is  ampli¬ 
fied  to  a  large  amplitude,  the  nonlinear  coupling  among  the 
frequency  components  of  the  wave  is  likely  to  occur  and 
should  be  addressed  in  the  future  work. 
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